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ABSTRACT 

We present multi-epoch observations of the radio nebula around the neutron star 
X-ray binary Circinus X-l made at 1.4 and 2.5 GHz with the Australia Telescope 
Compact Array between October 2000 and September 2004. The nebula can be seen 
as a result of the interaction between the jet from the system and the interstellar 
medium and it is likely that we are actually looking toward the central X-ray binary 
system through the jet-powered radio lobe. The study of the nebula thus offers a 
unique opportunity to estimate for the first time using calorimetry the energetics of 
a jet from an object clearly identified as a neutron star. An extensive discussion on 
the energetics of the complex is presented: a first approach is based on the minimum 
energy estimation, while a second one employs a self-similar model of the interaction 
between the jets and the surrounding medium. The results suggest an age for the 
nebula of ^ 10 5 years and a corresponding time-averaged jet power ^ 10 35 erg s — 1 . 
During periodic flaring episodes, the instantaneous jet power may reach values of 
similar magnitude to the X-ray luminosity. 

Key words: accretion, accretion discs - stars: individual: Circinus X-l - stars: vari- 
ables: other - ISM: jets and outflows - radio continuum: stars - radiation mechanisms: 
non-thermal - X-rays: stars. 



1 INTRODUCTION 

Circinus X-l is a very unusual X -ray binary system discov- 
ered in 1971 <|M argon et al.|l97lh . It und ergoes outbursts at 
X-ray jKaluzienski et alJll97rJ)~infrared JGlasslll978t 1199411 
and radio jwh^?an^^n7H977t lHavnes etHnTfl97'8T) wave- 
lengths with a period of 16.6 days, a fact that was inter- 
preted as enh anced accretion near the periastron of an el- 
liptic al orbit jMurdin et al.l Il98d: INicolson. Glass fc Feastl 
1980). Millisecond variability in X-ray led initially t o the 
classi fication of Cir X-l as a black hole candidate llToorl 
Il977f) . but discovery of type I X-ray bursts f rom the field 
of Cir X-l llTennant. Fabian fc Shafeil Il986allbl suggested 
identification of the compact object with a neutron star. 
The quasi-periodic oscillations (QPOs) and X-ray colors of 
the object on occasion exhibit behaviour similar to that of 
Z sources JShirev et al.l Il999h but at other times Cir X- 
1 resembles an atoll source ( Oosterbroek et all Il995f) and 



the source does not fit perfectly to either of these two 
flavou rs of low magnetic field neutron stars f see Ivan der Klisl 
( 1995)). Optical/IR obs ervations point to a low mass X -ray 
binary system (LMXB1 J Johnston. Fender fc Wull999l) . Re- 
cently, iBoutloukos. Wiinands fc van der Klisl ^OOeF found 
twin kHz QPOs in the power spectra offering more evi- 
dence for associating the compact object in Cir X-l with 
a neutron star. The orbit of the bi nary is not known 
but probably has a high eccentricity jMurdin et al.l ll98(J : 
INicolson et alJll98d ; Ffauris et al.lll999l). C ir X-l lies within 
a synchrotron nebula ljHOTn^s^taLl^ii^) and has variable 
radio flux densities at cm wavelengths. W hile the radio 
flares reached up to 1 Jy in the late 1970s iWhelan et alJ 
Il977t lHavnes et aljll978Ti . they have only been observed a t 



the tens of mjy level ever s i nce (e .g. IStewart et al.l il99ll) : 
iFender. Tzioumis fc Tudosd i2005t0 . The age of the system 
is unknown, almost all of the estimates being made within 
the framework of the "runaway binary" scenario that arose 
from the hypothetical associati on with the nearby supernova 
remnant (SNR) G321. 9-0.3 JClark. Parkinson fc Caswell 
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Il975l) . However, more recent Hubble Space Telescope (HST) 
observations have revealed no evidence for proper motion 
of Cir X-l, almost certainly ruling out such a relation 
jMignani et al.ll2002l) . leaving the age of Cir X-l extremely 
poorly constrained. Another important parameter, the dis- 
tance, has a very large range of possible values, rangin g from 
4 to 9 kpc JGoss fc Meboldll977t llaria et aljl2005h . This 
system also seems to harbo ur the most relativ istic outflow 
observed within our galaxy l|Fender et alJl2004h . The arcsec 
scale jet has a minimum bulk Lorentz factor of around 10 
and an inclination with respect to the line of sight of less 
than ~ 5 degrees, depending on the accepted distance. As a 
result, Cir X-l resembles, at least from a geometrical point 
of view, a down-scaled version of a BL Lac object. 



2 OBSERVATIONS 

We have observed Cir X-l over five epochs between October 
2000 and September 2004 at 1.4 and 2.5 GHz using the Aus- 
tralia Telescope Compact Array (ATCA). For better uv cov- 
erage, four different antenna configurations were employed, 
with baselines ranging from 31 up to 6000 m. Table 1 con- 
tains the observational log. The absolute flux density cali- 
bration was scaled with respect to PKS J 1939-6342 (PKS 
B1934-638). For phase calibration the nearby PMN J1524- 
5903 (B1520-58) was used. At each of the two frequencies 
the data were calib rated for the five epochs individ ually with 
MIRIAD software JSault. Teuben fc Wrighll995l) and then 
combined together into an image that covered this way the 
widest range of baselines permitted by ATCA. In our study 
we also made use of the publicly available 1 0.8 GHz Molonglo 
Observatory Synthesis Telescope (MOST) data (FITS image 
from April 3, 1992) which is part o f the Sydney Univer- 
sity M olonglo Sky Survey (SUMSS) dBock. Large fc Sadlerl 

The nearby projection on the sky of SNR G321.9-0.3 
(see Fig. 1) resulted in artifacts on the radio maps, degrad- 
ing the quality of the images (the ATCA full width at half 
power primary beam is 33 arcmin and 22 arcmin at 1.4 and 
2.5 GHz respectively). As a result of this and the diffuse na- 
ture of the Cir X-l nebula itself, the images and flux mea- 
surements of the nebula vary significantly with the chosen 
angular resolution (which we can select by adjusting the 
weighting on different baselines). Furthermore, for the to- 
tal flux measurements used in section 2.2 we weighted the 
short baselines most strongly in order to recover as much of 
the diffuse nebular emission as possible (while still trying to 
exclude artifacts from SNR G321.9-0.3). 



2.1 Radio maps 

In Fig. 1 we present the large-scale 0.8 GHz radio map of 
the Cir X-l region from SUMSS, indicating the nearby SNR 
G321.9-0.3. Fig. 2 is a montage of 0.8, 1.4 and 2.5 GHz maps 
of the nebula around Cir X-l, clearly showing structure and 
emission extending in the direction of the SNR, and finally 
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Figure 1. MOST total intensity 0.8 GHz map of the SNR G321.9- 
0.3 and radio nebula around Cir X-l. The contours are at -1, 1, 
1.4, 2, 2.8, 4, 5.6, 8, 11, 16, 23, 32, 45, 64 and 90 times the rms 
noise of 8.7 mjy/beam. 



Fig. 3 presents our most detailed 1.4 GHz image of the neb- 
ula, revealing considerable structure within the arcmin-scale 
jets. 

An arcmin jet-like str ucture within the neb ula (previ- 
ously observed in radio bv lStewart et all (^)93)) is clearly 
seen at 0.8, 1.4 and 2.5 GHz in Figs. 2 & 3. The jet exhibits 
asymmetry: the south-eastern side (with respect to the radio 
core) is slightly brighter and longer than the north-western 
one. 

As already noted, the arcsec scale jet of Cir X-l may 
be the most relativistic yet discovered in the galaxy, with 
an or ientation very close to the line of sight jFender et alJ 
2004) and a pos ition angle on the sky coincident with the 
arcmin scale jet iFender et aljll998Tl . It therefore seems rea- 
sonable to assume that the geometry of the jet is preserved 
at arcmin scales, so the observed asymmetry at 1.4 GHz can 
be explained in terms of an approaching (south-east) and a 
receding jet (north-west) affected by projection effects. The 
bends of the jet may be consequences of the precession af- 
fecting the disk of the X-ray binary to which the jets are very 
likely co upled. Such pr e cessio n of the jets is observed in SS 
433 Ce.g. lMargon et alJ J1979D 1. an X -ray binary harbouring 
relativistic jets interacting with a surrounding radio nebula, 
an object resembling Cir X-l. Interactions with clouds of 
higher density could be an alternative explanation for the 
change in direction of the jet an d this would appear to be 
observed at arcsec scales at least jFender et alJl2004|) . 

The expected ratio between the fluxes of the approach- 
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Table 1. Observational log. 



Epoch 


Time [h] 


Array configuration* 1 


Frequency [GHz] 


October 1, 2000 


9 


6A 


1.4, 


2.5 


October 14, 2000 


11 


6A 


1.4, 


2.5 


October 25, 2000 


12 


6C 


1.4, 


2.5 


August 3, 2001 


11 


1.5A 


1.4, 


2.5 


September 3, 2004 


10 


EW214 


1.4, 


2.4 



a A list of standard ATCA array configurations at: 

|http://www.narrabri.atnf.csiro.au/operations/array_cormgurations/configurations.html| 
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Figure 2. Left: MOST total intensity 0.8 GHz map of the radio nebula around Cir X-l. The contours are at -0.5, 0.5, 0.7, 1, 1.4, 2, 2.8, 
4, 5.6, 8, 11, 16, 23, 32, 45, 64 and 90 times the rms noise of 8.7 mjy/beam. The beam size is 51 X 43 arcsec 2 at PA=0°.0. Center: 
ATCA total intensity 1.4 GHz map of the radio nebula around Cir X-l. The contours are at -1, 1, 1.4, 2, 2.8, 4, 5.6, 8, 11, 16, 23, 28, 32, 
37, 45, 64 and 90 times the rms noise of 0.7 mjy/beam. The beam size is 26.2 X 23.5 arcsec 2 at PA=4°.9. Larger than nominal ATCA 
beam size was used in order to show maximum of details at all scales. Right: ATCA total intensity 2.5 GHz map of the radio nebula 
around Cir X-l. The contours are at -0.7, 0.7, 1, 1.4, 2, 2.8, 4, 5.6, 8, 11, 16, 23, 32, 45, 64 and 90 times the rms noise of 6.1 mjy/beam. 
The beam size is 33.7 X 29.4 arcsec 2 at PA=-6°.8. Larger than nominal ATCA beam size was used in order to show maximum of details 
at all scales. 



ing and receding continuous jets, for an inclination with re- 
spect to the line of sight of 5-15 degrees and a velocity ratio 
f3 = v/c = 0.05, 0.10, 0.15) is (1.3, 1.6, 2.1). An image plane 
fit to the central region of the Cir X-l complex in the radio 
map in Fig. 3 (the inner 40 x 40 arcsec 2 ) shows that the sim- 
plest appropriate model is composed of an unresolved source 
coinciding on the sky with the X-ray binary (at the position 
given by Fender et al. 1998) and an extended source centred 
some 2 arcsec to the south-east. If we assume that the SE- 
NW asymmetry is due to Doppler boosting, it follows that 
the radio map at 1.4 GHz is consistent with a large-scale 
jet with a f3 factor of at least 0.1. The maps at the other 
frequencies have a poorer resolution and cannot add useful 
supplementary information in this respect. As a caveat we 
should note that it is clear from the results presented in 
Fender et al. (2004) that complex and varying interactions 
do take place within the inner few arcsec, and there may be 
a significant contribution to the SE-NW asymmetry from 
e.g. the environment into which the jet is propagating. 



An interesting feature present at all three frequencies 
is "knot A" in the north-east (Figs. 2 & 3). Its nature is 
unknown. It could be a region of higher density that was 
energised by the passage of a shock originating in the center 
of the complex. The blob is not detected at higher frequen- 
cies (October, 2000 ATCA 4.8 and 8.6 GHz) therefore an 
identification with a background radio source can probably 
be ruled out, since the lack of detection at these frequencies 
is likely due to a diffuse structure being resolved out. 

An enhanced radio emitting region, "knot B" , is evident 
at 0.8 and 1.4 GHz (Figs. 2 & 3). This could be interpreted 
as the interaction site where the approaching jet impinges on 
the interstellar medium. However, projection effects might 
be at work and could complicate such a simple explanation. 

The "knot C" region (Fig. 2) also appears at all three 
frequencies discussed here and looks rather isolated at 0.8 
GHz and 1.4 GHz. The nature of this is unknown. Its pres- 
ence was easier to accommodate in the "runaway binary" 
scenario, but the HST observations (see section 1) reopened 
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Figure 3. ATCA total intensity 1.4 GHz map of the central part of radio nebula around Cir X-l with data from August 2001 only 
(i.e. the data set doesn't contain the shorter baselines from epoch September 3, 2004 and is therefore less sensitive to the diffuse radio 
emission, neither the data from October 2000 which have a poorer u-v coverage due to frequency switching). The contours are at -1, 1, 
1.4, 2, 2.8, 4, 5.6, 8, 11, 16, 23, 32, 45, 64 and 90 times the rms noise of 2.0 mjy/beam. The beam size is 18.1 X 16.5 arcsec 2 at PA=-7°.l. 
Larger than nominal ATCA beam size was used in order to show maximum of details at all scales. 



the issue of its origin. It may be a relic of a previous epoch 
of strong jet activity. 

The likely general picture of the radio nebula emerging 
is that of a system in which the energy and matter are sup- 
plied from the central source via the jets. In this respect, 
we propose that the radio nebula is formed at the interac- 
tion site between the jet, inclined at a low angle to the line 
of sight, and the interstellar medium (Fig. 4). Therefore, 
uniquely in the case of this system, we are observing the 
core of the binary through one of the radio lobes. 

2.2 Spectral index 

Even though the radio monitoring of Cir X-l was sparse in 
the last three decades, apparently there is a secular decrease 
in the radio flux levels measured in the active as well as in 
the quiescent states, wi t h short periods of enhanced activ- 
ity JWhelan et alJll977t lHavnes et al"lll978l : iNicolson etafl 




radio lobes 



Figure 4. The proposed geometry of the Cir X-l complex (not 
to scale) based on analogy with the W50/SS433 complex. 



198(1 JPreston et alJll98A IStewart et alJll99ll Il993t iFenderl 
1997tlFender et alJl998ll2005l) . The ASM /Rossi X-ray Tim- 
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Figure 5. The overall spectrum of Cir X-l complex in the qui- 
escent state. 



ing Explorer (RXTE) 2 data from the last ~ 10 yr show a 
similar behaviour in X-ray band. Fig. 5 presents the quies- 
cent overall spectrum of the Cir X-l complex (i.e. binary 
system plus radio nebula) from combined ATCA and MOST 
data (see section 2). This is the first low-frequ ency spectrum 
reported for this object since the lat e '70s iWhelan et alJ 
Il977l: lHavnes et aljll973) . early '80s (lHavnes et al.lll98fj) . 
The spectral index obtained here, a ~ —0.5 (F v oc v a ), is 
within the errors identical to earlier values, with the flux 
level scaled down by ~ 300 mjy at 1.4 GHz in comparison 
with the measurements (notably at a sligh tly different res- 
olution) made in 1977 (lHavnes et alJll978l) . This is almost 
certainly due to the secular decline of the radio core and we 
find no evidence for a change in the radio luminosity of the 
nebula. 

The spectral index map of the radio nebula, made with 
the combined ATCA data at 1.4 and 2.5 GHz, is presented 
in Fig. 6. A larger than nominal beam was used in order to 
increase the signal to noise ratio of the extended emission 
(nevertheless, due to the simultaneity of the observations at 
the two frequencies, we need to bear in mind that the u-v 
coverage is not identical). Maps made with better resolution, 
which at least in principle should be qualitatively correct, 
are compatible with the features seen in Fig. 6. Namely, there 
is a tendency towards a flatter spectrum in the north-eastern 
part of the nebula in a zone coincident with an enhanced ra- 
dio emission region ("knot A" in Figs. 2 & 3). This could be 
a site of augmented electron acceleration. Also, at a lower 
confidence level, the spectrum seems to flatten at the edges 
of the nebula. If real, this could be interpreted in terms of 
free-free emission or synchrotron radiation from enhanced 
accelerated particles at the interaction site with the inter- 
stellar medium. The "knot C" (Fig. 2) doesn't appear on 
the map due to our chosen cut-off in the confidence level of 
the features represented. 




15°31 m 10 s 21 m s 50' 40 s 30 s 20 s 10 = 
RA (J2000) 

Figure 6. The spectral index map of the radio nebula around 
Cir X-l between 1.4 and 2.5 GHz. The beam size is 139 X 108 
arcsec 2 at PA=-22° .4 



3 ENERGETICS 
3.1 Minimum energy 

Assuming a source of volume V, specific luminosity L v , spec- 
trum of the form L v oc v a and considering the radiation 
to be of synchrotron origin produced by a particle pop- 
ulation with the energy spectrum N(E)dE — NoE~ p dE, 
it is possible to estimate the minimum energy associated 
with the source. The spectral index a is related to the par- 
ticle distribution index p by p — 1 — 2a, N(E)dE is the 
number density of particles in the interval [E, E + dE] and 
No is a proportionality constant. Two main pathways have 
been followed in literature for determining the relevant for- 
mulas: one usin g a fixed interval in frequency in t he radi- 
ation spectrum feurbidgd Il956t iPacholcz vkl 1 1 9 701 : iLongair 
Il994l ), the other in the particle energy spectrumfe.g. [ Pohll 

BlirJ 



1 19931: iBrunetti. Setti fc Comastrilll997t iPfrommer fc E nfili 
2004). We will adopt the former method and closely follow 
lLongairl (ll994l). 



The internal energy associated with the observed syn- 
chrotron emission can be split in two components: E e the 
energy of particles (i.e. electrons and protons) and Eb the 
energy of fields (i.e. magnetic): 



E c = C e (l + V)L V B 
E B = C B fVB 2 , 



-3/2 



(1) 
(2) 



http: / /xte. mit.edu/ 



where B is the magnetic flux density, rj is the ratio between 
the energies in protons and electrons, / is a parameter tak- 
ing into account the uncertainties in the estimation of the 
volume V of the emitting region (it is not what usually is 
referred to as the volume filling factor because we allowed 
it to be higher than 1, covering this way also the situations 
when an underestimation of the volume is made; this is very 
relevant in the case of Cir X-l because there may be large 
projection effects along the line of sight). Cb is a numerical 
constant related to the magnetic constant /10 , while C e in- 
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corporates the dependences on t he radiation sp ectrum and 
particle energy distribution (e.g. iLongairl ]l994h '). 

The sum of energies in eqs. (1) and (2) has a minimum 
with respect to the magnetic flux density at : 



D,. <m = 2.:; I 10 -■ C '(I -ri) 2/7 .r V7 (^~~ : 



VWHz- 1 / 



2 st2/7( 
e 

2/7 



T, 



(3) 



thus the minimum energy necessary to explain the observed 
synchrotron emission reads: 

V \ 3/7 



E 



= 4.97 x 1(T C, 



t /r {1 + vr/ r f/ r(V_y 



X (i 



,WHz- 



4/7 



(4) 



-^min 

and E m i n correspond to a particular state of the system 
in which the energy in magnetic field is three quarters of the 
energy in particles. 

Averaging the energy loss rate by synchrotron radiation 
over an isotropic distribution of pitch angles, the lifetime of 
an electron emitting at the peak of the synchrotron spectrum 
is: 



t 



8.5 x 10 s 



HzJ I T j 



-3/2 



(5) 



3.2 Application to Cir X-l 

In order to estimate the minimum energy requirements for 
Cir X-l we use eq. (4). The n ebula was detected between 
0.4 GHz dWhelan et al] lll977li : considering that at low fre- 
quencies the radiation from the diffuse nebula dominates) 
and 2.5 GHz (this work), therefore we make the assump- 
tion that within this frequency interval the spectrum is of 
synchrotron origin and has the slope derived in section 2.2, 
a = —0.5. The energy E e associated with the particles is 
dominated by the lowest energy particles and so the above 
upper limit in the frequency domain is not essential in the 
calculation. The volume of the emitting region can be very 
roughly estimated in this particular case by assuming the 
radio lobe has a sp herical geometry: for a 250 arcsec radius 
object at 4.1 kpc Jlaria et al] I2005T) the resulting volume 
(both radio lobes taken into account) is 2.2 x 10 52 m 3 . We 
further assume that the synchrotron radiative processes are 
dominated by electrons and therefore the ratio between the 
energy in protons and electrons r\ = 0. Applying eq. (4), 
the minimum energy (in the comoving frame) required to 
account for the observed radio emission, under the assump- 
tions made above is 8.1 x 10 46 erg. This value is of the same 
order of magni tude as the one infe r red for the W50 nebul a 
around SS433 JPubner et al.lll998l: iMoldowan et al.ll2005l) . 
This raises the possibility that Cir X-l might not be a jet- 
powered radio nebula o nly, but mayb e a distorted SNR as 
suggested for W50 (e.g. iGreenl ]2004fl 'l. The corresponding 
magnetic field derived from eq. (3) is 6.3 uG. The electron 
lifetime (eq. (5)) at 2.5 GHz is then 3.4 x 10 7 yr. 

This estimate does not include energy stored in bulk 
motions or in heating of protons and so could be a significant 
underestimate. Errors in the evaluation of the volume (/ g 
[10 -2 , 10 2 ]) would modify the minimum energy by not more 
than a factor 10. 
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Figure 7. The averaged jet power of the radio nebula around 
Cir X-l as a function of its age and the duty cycle of the binary 
system. 



Assuming the energy in the nebula has not been in- 
jected continuously (i.e. along the full orbit of the X-ray bi- 
nary), Fig. 7 shows the dependence of the averaged jet power 
as a function of the age of nebula and the duty cycle of 



X-ray 


flares correlate with the radio flares J 


Whclan 


et al. 


Il977t 


Havnes et alJ 


1978:lThomas et all 


1978; 


lFenderlll997l) 



suggesting enhanced accretion onto the compact object and 
ejection of matter via the jets. It is realistic to consider that 
the injection of energy from the binary to the radio nebula 
is occurring mainly during the flare states and to identify 
the duty cycle in Fig. 7 with the duration of the flares. The 
solid (lowest) line in Fig. 7 corresponds to a duty cycle of 
unity. For an active jet phase of less than one day (corre- 
sponding to a duty cycle 0.06), and an age of the nebula 
of 10 4 — 10 s yr , the minimum jet power would be of the 
order of 10 36 - 10 37 erg s _1 . 

Fig. 8 presents the ASM/RXTE light curve of the cen- 
tral X-ray source of the Cir X-l complex. Over the period 
of the X-ray data plotted, the averaged X-ray power (for a 
distance of 4.1 kpc) is 4.3 x 10 37 erg s _1 . Such a high power 
is typical for most of the Z sources, which are accreting at 
a large fraction of the Eddington limit. 

Making the simple assumption that in the flare state 
the X-ray light curve has a square shape and the increase 
in amplitude from the quiescent level is of 50 ASM counts 
s _1 (a more or less typical value in the last few years), then 
the additional power in X-rays during the flare is of the 
order 3.4 x 10 37 erg s _1 . In any case, these estimates suggest 
that the jet power during flaring phases may be a significant 
fraction of the X-ray power. 



4 AGE OF THE RADIO NEBULA 

Constraining the age of the nebula is a very difficult task 
given the fact that so little is known about the Cir X-l 
complex. Only very weak constraints can be imposed. 

The size of the nebula can offer a rough estimation of the 
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Figure 8. 2-10 kev ASM/RXTE X-ray light curve of Cir X-l 
between February 1996 - July 2005. 



age. Oversimplifying the problem and considering a constant 
expansion of the nebula with an average velocity in the plane 
of the sky v, the time required to achieve its present size 
(assuming a distance of 4.1 kpc) is: 



''expansion 



= 2.0 x 10* 



km ; 



yr, 



(6) 



where 6 is the angular size. For v = 
(10; 50; 100; 500; 1000 km s -1 ), the expansion time is 
tcxpansion = (4 x 10 5 ;8 x 10 4 ;4 x 10 4 ;8 x 10 3 ;4 x 10 3 yr). 
Therefore if the nebula is expanding with a velocity 
comparab le to that of the jet powered nebula around Cyg 
X-l JCallo et al]l2005ft . that is tens up to a few hundreds 
km s _1 , its age is likely to be less than ~ 10 5 yr. 



5 A SELF-SIMILAR FLUID MODEL 

In order to attempt to understand the energetics and evolu- 
tion of the nebula better, we have developed a simple self- 
similar fluid model for the jet-powered nebula. 

We identify the radio nebula of Cir X-l with the ra- 
dio synchrotron lobes inflated by the jets. This situation is 
analogous to the model for the lobes of radio galaxies de- 
scribed bv lKaiser fc Alexanded dl997t) . In the following we 
develop a simplified version of this model for Cir X-l. The 
main assumptions of the model are: the jets are in pressure 
equilibrium with the lobes they inflate. The energy trans- 
port rate of a single jet or "jet power", Qo, is constant over 
the lifetime of the jet. The jet may go through phases of out- 
bursts and quiescence as long as the duration of these phases 
is short compared with the overall jet lifetime. In this case, 
Qo is the jet power averaged over the lifetime of the jet. The 
expansion of the lobe is confined by the ram-pressure of the 
rece ding external gas in all dir ections. Under these condi- 
tions |Kais£r_&^Alex£mde]] 1119971) showed that the expansion 
of the lobe is self-similar, e.g. the ratio R = Lj (2r), where 
L is the length of the lobe in the direction of jet propaga- 
tion and r is the radius of the lobe perpendicular to the jet 
measured at a fixed fraction of L, say L/2, is a constant. 
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Figure 9. Cylindrical radio lobe seen from the side. In Cir X-l 
our line of sight is very close to the jet axis. The relative size of 
the head region is exaggerated. 



We now simplify the model by assuming a cylindrical 
geometry for the lobe with the axis of the cylinder aligned 
with the jet. Figure [5] shows a schematic of the lobe ge- 
ometry. The jet ends in a high- pressure region, the "head" , 
where the energy transported by the jet is thermalized. The 
pressure in the head, ph, mainly drives the forward expan- 
sion of the lobe, which is confined by the ram-pressure of the 
external medium of mass density p. This can be expressed 



Ph 



L p, 



(7) 



where we assume that the density of the external medium, 
p, is constant. The sideways expansion of the lobe is mainly 
driven by the pressure in the lobe, p\, and so 



Pi 



• 2 

r p. 



The ratio of the two pressures is then: 

2 



Ph 
Pi 



4ir 



(8) 



(9) 



where we have used the condition of self-similar expansion, 
i.e. r — Lj (2R) with R = constant. 

The volume of the head region, Vh, will usually be small 
compared to the volume of lobe, V\ and so Vh/Vi = Rv <C 1 
jKaiser fc Alexanderl ll999T) . The self-similar expansion of 
the lobe suggests that Rv = constant, but we will see in 
the following that as long as Rv € 1/ (4^? 2 ) anv time de- 
pendence of Rv is unimportant. 

The total internal energy of the lobe, U, contained in 
V\ changes according to: 



dU = 



■ (Mdpi + pidVi) = Qo&t - pidVi - phdVi, (10) 



7-1 

where 7 is the adiabatic index of the lobe material. The last 
term describes the expansion work of the head region. Using 
the results of the discussion above we find: 



p h dU h = 4R 2 R v pidVi. 



(11) 



Thus we can eliminate all quantities referring to the head 
region from equation ill 01 . From equation ||§J and remem- 
bering that r — Lj (2R) we get: 



dt = -LJ£dL. 

2RV Pl 



(12) 



Finally, the self-similar expansion of the lobe and Rv = 
constant imply 

L ^ 3 



where Lo is an arbitrary scale length 



(13) 
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Substituting equations [11 1 1 , I12H and 11311 into equation 
11011 yields after re-arranging 



-3 [7 + (7 - 1) 4R 2 R V ] S.(14) 



The solution of this equation is: 
L \ " 4/3 



Pi = Po 
with 

po = 



(15) 



7-1 



QoLq 



2-R{3[ 7 +(7-l)4R 2 i?v]-4/3} V 



2/3 



\/l> f -(16) 



From the observations we can determine r and an esti- 
mate for p\. We then have to assume a value for the aspect 
ratio R, where 1 ^ R ^ 5 is reasonable. With this we can de- 
termine the length of the lobe L. Since we are free to choose 
Lq, it is convenient to set Lo = L. Vq then follows from Lq 
and R and an assumption about Rv ■ If Rv < 1/ (4i? 2 ), 
then Vo ~ nr 2 Lo = kLq/ (4i? 2 ) and we can neglect the 
term involving Rv in the eq.(21). With a typical value for p 
appropriate for the external medium we can then determine 
Qo, the time-averaged jet power. 

Finally, the age of the jet/lobe system can be found 
from equation (1121 and is, for our choice Lq = L, 



as a galactic analogue of the radio lobes in AGNs: the jet 
interacts with the interstellar medium and creates the syn- 
chrotron radio emission. Uniquely in this case, we are prob- 
ably seeing the central source through the jet-inflated radio 
lobe, as the jet appears to be aligned very close to the line of 
sight. In this respect, the Cir X-l nebula may be an "end-on" 
analogue of the SS 433/W50 com plex. The behavior of the 
small scale jet dFender et al ] |2004l) suggests that the energy 
that powers up the nebula is supplied intermittently from 
the central binary system via the jets. Even though hard 
to strongly constrain in value, the amount of energy trans- 
ferred to the nebula is important and should be taken into 
account whenever an energy balance of the whole system is 
required. The calculations suggest an averaged jet power of 
at least 10 35 erg s~ x for a preferred age of the nebula of ^ 
10 years and it may be that the jet power becomes compa- 
rable to the X-ray luminosity during outbursts. Transfer of 
considerable energy via jets is not rare, th e case of SS433 be- 
ing well studied (e.g. lDubner et aD l|l998f0 : moreover such a 
process might be more co mmon than previou sly suspected, 
as pointed out recently bv lGallo et al] i2005l) . 

Finally, we note that this is one of the best estimates 
to date for the power in jets from a neutron star system. 
Such measurements may be key in attempts to quantify the 
disc-jet coupling not only in ot her neutron star, but also in 
black hole systems of all scales llKording et al.ll2006l) . 



t ■■ 



(17) 



We assume the lobe is populated by relativistic elec- 
trons and therefore the adiabatic index 7 = 4/3. pn can be 
estimated from the minimum energy requirements as dis- 
cussed in section 3 and using eq. (14) a minimum lobe pres- 
sure pi is obtained. For an angular size of the lobe of 250 
arcsec at 4.1 kpc, R=4, Rv = 0.01 and a number density of 
the external medium of 2 cm -3 , the resulting long-term av- 
eraged jet power is 4x 10 35 erg s _1 while the age is 2x 10 4 yr. 
For lower duty cycles, as is very probably the case for Cir 
X-l, the jet power estimated here could be up to two or- 
ders of magnitude higher. The orders of magnitude of the 
jet power and age of the nebula are weakly dependent (i.e. 
do not affect the statements we make) on different reason- 
able values assumed for the parameters of the model. For 
instance, for a distance d a factor two higher than the one 
we preferred, the jet power becomes 7.5 x 10 34 erg s _1 and 
the age 1.1 x 10 5 yr. 

With an averaged jet power of ~ 10 35 erg s _1 and an 
age of the nebula of less than ~ 10 yr, these results are 
compatible with the ones presented in section 3 (see Fig. 7). 
Clearly, it is not unphysical that the jet could have inflated 
the nebula. 



6 CONCLUSIONS 

We presented ATCA radio maps of the nebula around Cir 
X-l at 1.4 and 2.5 GHz. Combining them with publicly avail- 
able MOST data we obtained the quiescent spectrum for the 
complex. The slope of the spectrum is within the error iden- 
tical to the ones obtained in the 1970's-1980's but is scaled 
down by around 300 mjy at 1.4 GHz, compatible with the 
observed secular evolution of the central source, both in ra- 
dio and X-rays. As a whole, the radio nebula can be seen 
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